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by 
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SUMMARY 

A simple model of the interplanetary magnetic field is described 
and solved analytically. In this model space is divided into three 
regions by two concentric spheres. Conductivities (with one exception) 
are assumed to be isotropic and constant in each region, and flow 
velocities are regular and prescribed. The innermost region rotates 
rigidly around its center, the intermediate region contains a com- 
pressible fluid flowing radially outward at a constant velocity (an 
idealization of the solar wind), and the outer region is at rest. The 
magnetic field originates at point sources at the origin and possibly in a 
uniform field at infinity. With these assumptions methods are described 
for finding the field in the general case, and also in the limit when all 
conductivities are very high. As an example the case in which the field's 
source is a point dipole aligned with the axis of rotation is solved in 
some detail. Part II of this study, published as a separate Technical 
Note, considers the cosmic ray anisotropy. 
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A SIMPLE MODEL OF THE 
INTERPLANETARY MAGNETIC FIELD, 

PART I: CALCULATION OF THE MAGNETIC FIELD* 

by 

David Sternt 

Goddard Space Flight Center 


INTRODUCTION 

Experimental evidence indicates that an interplanetary magnetic field originates at the sun 
and extends at least to the earth T s orbit, possibly much further. Many features of this field are 
still uncertain, but two of its main properties have been predicted theoretically and seem to agree 
with observation. Both may be regarded as manifestations of the fact that a highly conducting 
fluid — here the solar wind emanating radially from the sun — tends to impart its motion to lines 
of force embedded in it. It was theorized that the solar wind stretched the lines of force, render- 
ing them almost radial and causing the field intensity B to fall off less rapidly than it would other- 
wise (see for example Reference 1). In addition to this the field was expected to be twisted by 
solar rotation into an Archimedean spiral. This point was noted first by Chapman (Reference 2) 
who observed that the locus of a particle stream constantly emitted from a point on the sun is, at 
any time, such a spiral. (Because the same locus is described by droplets from a rotating 
sprinkler, this is sometimes called the "garden-hose effect.") A line of force drawn out by a 
stream of particles would also follow such a spiral, and it was argued that similar twisting occurs 
in any field originating in the rotating sun. Parker gave a formal proof of this, assuming that the 
magnetic field is parallel to the velocity field as seen from a frame of reference co-rotating with 
the sun (Reference 3). The effect has also been deduced from experimental data, from the arrival 
direction of solar flare particles (Reference 4), and from direct observation (Reference 5) by 
Mariner 11(1962 apl). For these cases the "garden-hose angle" between B and the radial direction 
from the sun was of the order of 45 degrees. 

In this work a simple model of the interplanetary magnetic field will be investigated, first in 
the limiting case of a perfectly conducting fluid and then for finite, isotropic, and homogeneous 
conductivity. In Part II of this study the cosmic ray anisotropy is investigated (Reference 6). The 
description of the model follows. Space is assumed to be divided into three regions by two 


♦This report supersedes Goddard Space Flight Center document X-640-63“l62, August 1963. 

^The major part of this work was performed while the author held a National Academy of Sciences — National Research Council Post — 
Doctoral Resident Research Associateship. 
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Figure 1— The division of space into three regions. 


concentric spheres of radii R Q and Rj (Figure 1). 
Region I, the innermost region, is assumed to 
rotate rigidly with angular velocity <u. This 
region contains the source of the magnetic field, 
which will be assumed to be concentrated at the 
origin. Because of the tedious calculations 
involved, only the case of a dipole source with 
moment M will be treated in any detail. If this 
dipole is inclined at an angle k to the axis of 
solar rotation, its axial component may be re- 
garded as an idealization of the main solar 
field, and its equatorial component as that of an 
active region co-rotating with the sun. Region 
II, between the spheres, contains a compressible 
conducting fluid flowing radially outward with 
constant velocity u. Finally, in region III, 
which extends to infinity, no motion takes place. 
Region I actually represents the sun, and region 
II the space swept by the solar wind. The 
equation connecting the magnetic field B, the 
electric field E, the velocity v, and the conduc- 
tivity C7 iS 


curl B - fi Q ct[E + (v x B)] . (1) 


We shall be interested in stationary solutions, and particularly in the case when a is large and 
>:> *o * 

It should be borne in mind that the preceding is a gross simplification of the actual situation. 

To stress this point, everything approximated and neglected will now be listed: 

1. It has not been established that the solar dipole plays a major role in creating the inter- 
planetary field. Certainly, the source of the field is quite complex. 

2. The solar wind is not an ordinary conducting fluid but a nearly collisionless plasma, con- 
ducting very well along the magnetic field but much less across it. Unfortunately, since 
the direction of the conduction anisotropy depends on the magnetic field, taking it into 
account makes the equation of conduction nonlinear. Except for one case, therefore, the 
conductivity will be assumed to be isotropic. 

3. The flow and.field assumed here are laminar and regular, although observation indicates 
a large irregular component. The model developed here thus represents only the effects 
of the average interplanetary field and does not include the turbulent component. 
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4. The assumed sharp boundaries are only a convenient approximation to the actual ones, 
which are not too well known at the present time. 

5. In any problem of this sort, the velocity v generally is not determined a priori but is 
solved simultaneously with B, by using the hydromagnetic flow equation (see for example 
Reference 7). This equation is nonlinear and with this approach the term vxB in Equation 1 
becomes nonlinear too. In the vicinity of the earth, of course, the mass flow dictates the 
magnetic field because of its much higher energy density; nevertheless, in the vicinity of 
what corresponds to the outer sphere of this model, the flow may be considerably distorted 
by the field. 

Unfortunately, a more realistic model would be very hard to solve analytically. It is hoped, 
however, that the results obtained here will give some qualitative insight into the behavior of the 
actual interplanetary field. 

INFINITE CONDUCTIVITY 

The field produced when has been derived by Parker (Reference 3), with a rotating 
frame of reference. It will be derived here in a somewhat more conventional way. In general, if 
the conductivity tends to infinity and curl B does not, 

E = - (v x B) . (2) 

Taking the curl in region II gives, in spherical coordinates (r, < 9 , <£), 

<?B 

dt = ucurl («* • (3) 


Utilizing 


gives 


div B - 0 


(4) 


where X stands for either r 2 B r 
and 


dX 
dt “ 


rB 0 , Or rB d . 


dX _ 

u Jr ~ 0 - (5) 

Under rotational symmetry x is independent of time 


r ~ 2 , 

( 6 ) 

US )*' 1 . 

m 
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and, by using Equation 4, 


B* = 0 . (8) 

E d is continuous on the surface r = R 0 . Just inside the boundary, by Equation 2, 

E e = “ ^R 0 sin 0 B r , 


and because B r is continuous its value just outside the boundary may be used. The continuity of 
E e then gives 




( 9 ) 


so that the tangent of the "garden-hose angle" i 

tan X " IT 

In general the source of the field rotates: 

B( r , 0 , 0, t ) 

K 

d 

at 

Then Equation 5 becomes 

ax _u 

aX* + co 


■■ -(^)rsin 0. ( 10 ) 

= B( r , B, A.) , (H) 

= 4> - cot , ( 12 ) 

= ~ ' (13) 

7=0, (14) 


and its solution is 


X = 



- cot 


COT \ 

“u~ ) • 


(15) 


Thus 


B r - r ' 2 (S) expim (<£-&*- ^ 7 ) , 


(16) 
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and the continuity of E e when r = r 0 gives 




(ll?) Sin 6 ( 6 ) eX P im {$-<**- ir) • 


(17) 


For a more detailed solution it is generally better to solve Equation 1 for finite conductivity and 
then investigate the behavior of the solution when a gets large. 


FINITE CONDUCTIVITY 

We will now turn to solving the problem for an arbitrary o-. It will be assumed that <y is 
uniform in each region and takes the values a 9 a 2 , and a 3 in regions I, II, and III, respectively. 
The following theorem is found useful: If a vector field B satisfies Equation 4, it may be uniquely 
resolved in the form 


B - curl (//j r + curl curl 0 2 r 


(18) 


Following Elsasser (Reference 8), we will term the component fields the toroidal and poloidal 
components, respectively. The theorem was proved first for rotational symmetry by Lust and 
Schliiter (Reference 9) and for the general case by Backus (Reference 10). In general, the follow- 
ing identities hold (Reference 11): 


curler = gradi p x r , 


curl curl i/r 


grad (<//r) - r V 2 \p 


By using Equation 18, a vector potential may be defined: 


A = \p x r + curl \p 2 r , 
B = curl A . 


A satisfies the gauge condition 


(19) 

( 20 ) 


( 21 ) 

( 22 ) 


divA = r-* 4~r (^, r 3 ) 


and the electric field E may be expressed as 


E 


. , ^ 
g rad '/ J 0 ■ dF • 


(23) 
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is 


The problem thus reduces to solving Equation 1 for three unknown scalars , \ p v and \p 2 . The 
derivation of a general solution tends to be tedious and therefore will only be outlined in this 
section. In each of the three regions Equation 1 can be put into the form 

grad^j + fq i + curlr£. = 0 ( 24 ) 

(i is the region's index). Taking the curl and applying the uniqueness of the resolution (Equa- 
tion 18) gives 


from which 


curl rr/ i - curl curl r £. - 0 , 


Vi ~ Vi (O , 
Cl = Ci(r) , 


and, after substitution in Equation 24, 


Si = Si (O • 


In general the expressions 77 . and £. will be functions of </q and \p 2 , which are also defined as 
arbitrary functions of r . If these functions can be chosen so that rj i vanishes, Equation 24 gives 


- const - C. 


Region I 

In region I 


v - (o x t - \ sin 6 . 

The following relations hold for v and any vector A : 


V 

- curl (oor 2 ) , 



curl v 

~ 2 CO , 



(A • V) v 

II 

e 

X 

> 



(v • V) A 

( dA r 

= <o x A + <ull r 

dk e 

+ , 'a 

+ l 4> d4> ) 


(25) 
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(for example see Reference 12, Equation 1.4.3). To resolve v x B in the manner of Equation 18 
note that 


v x B - v x curl A 

= grad (A ■ v ) - (Ax curl v ) - (A • V) v - (v • V) A 

t a A r «a*\ 

- grad (A • v) - &>li r -ftp +i $ d( p +i^ d( p \ , 

/ dp 2 \ / a^ 2 \ 

grad (cur sin 6 -ftf I - ru -jp - curl l r cu -jp I • 


B = 


(26) 


The rest of Equation 1 is easily resolved in the prescribed manner by using Equations 18 and 21- 
23, which give 


d 

Si O) = TF[Pi r ) + ^O 0- ! '/'o + Mo 0- ! " r sind ~W ’ 




(dpi 

dpi 

Vi ( r ) 

b 

0 

11 


+ OJ ~d$ 



( d P 2 

dp 2 


11 

"P 

0 

Ut 

+ “ ~Sp 


(27) 

(28) 
(29) 


The option to add an arbitrary radial function to<^, and p 2 i- s used to make both Vi (O and £1 (O 
vanish, so that 


d d ^2 

dTi'Pi r ) + Mo CT i Po + Mo 0- ! <ar sind ~m - C 1 ■ (30) 

So far no use has been made of the fact that the field's source co-rotates with region I. There- 
fore the equations hold even when the field originates, say, in region III as a fixed "interstellar 
field." If Equations 11-13 are introduced, 


V 2 Pi = 0 , 

(31) 

v*p 2 = 0 . 

<3 2) 


If the field's source is assumed to be concentrated at the origin, it is useful to expand p t andp 2 
in spherical harmonics; the expansion of P 2 then has a singularity at the origin corresponding to 
the source of the field. For instance, if this source is a dipole with moment M, inclined at an 
angle « to the rotation axis, then 


'Pi = a nm p n Pn” 1 (g) ex P im (<t>- cot) , 

nm 


( 33 ) 
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( 34 ) 


y * b nm p n P n m (O') exp im (<j6 - cot ) + [cos /< cos (9 + sin k sin 6 cos (4> ~ cot )] 


where 


r j Mo M 

p = -p- and /3 = — 

p R o p 47 tR 0 2 


Region II 

In this region 

v = ui r 

The contribution of the toroidal component to v x B is 

u(i r x curl i r r ) = u|j r x (gradi/^ r x i r jj 

= u jgrad (i//j r) - i r 'g~ ('/'j r)J 

and the contribution of the poloidal component is, by Equation 19, 

u i r X grad r)J = - curl r ^ 37 (i/> 2 r)J 

All the other terms of Equation 1 are the same as for region I. Therefore 


( 35 ) 


( 36 ) 


Si (O 

d , . 

- yPl r ) + M 0 a 2 'Po ~ u Mo a 2^1 r ’ 

( 37 ) 

v 2 (O 

U d .1 

- M 0 ct 2 pr + 7 Ti (Pi r )J - v Pi ■ 

( 38 ) 

s 2 (O 

= Mi^ 2 [-gr + TJti'Pi r )] - va ^a • 

( 39 ) 


As in region I, \p 1 and <f> 2 are chosen so that and S 2 vanish and g 2 equals a constant, c 2 . From 
Equations 11-13 the angular part of i/q is now expanded in spherical harmonics: 


'Pi = ^G n m ( r ) p " (0) exp im 


(40) 
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Inserting this into Equation 38 and using the independence of spherical harmonics gives, for any 
n and m , 


1 d 2 , 
r df 2 ( rG n 


■M 0 ct 2 


, . f n(n + 1) "1 

( rG nm ) - [ ~2 ima, Mo O- 2 Jl 


= 0 


(41) 


Multiplying by rR 0 2 and defining 


y nm(M) - P G wiP) 


gives 


U Mo r o y n m + [i Mo °-2 " R o 2 m " n (n + 1) p' 2 ] y nm - 0 


(42) 


The magnetic Reynolds number associated with radial outflow (in region II) may be defined by 


2a - M 0 °- 2 R o u . 


and that associated with rotation by 


2a 2o> = • 


Substituting 


y 


nm 


U nm eX P °'P 


and defining complex Reynolds numbers 2a 2c of order m, 

a 2c = " 2ima 2<» > 


we have 


“ % [ a 2 c +n(n + l)M‘ 2 ] = 0 . 


(43) 


With operators £ n defined by 


£ n = p-o> + « (p 3 


(44) 
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the general solution of Equation 43 may be written (Reference 13, p. 337, Equation 256) 


U nm - A 1 eX P a 2c P + A 2 K eX P " a 2c P 


= A l U n ,al(P) + A 2 (-!)” U ™2 <<°> ’ 

(45) 

By defining 


e . - p~ 1 u . exp ap , 

& nmi r nun ^ r 


with i = 1, 2, we can write the general forms of and ip 2 in region II as 


'Pi = 2^ P ” (0)eXpim( ^" a,t) [ a ™'l g ™”l ( ^ ) + a '>”2 g ™2 (/:,) ] * 

nm 

(46) 

'Pi = P n m eX P im & ~ ** ) [ b ™l gnml <-P) + b nm2 Snm2 (<°>] ' 

nm 

(47) 

Region III 


In this region Equation 1 has novxB term. Consequently 


d , . 

£3 ( r > " TF ('Pi r ) + Po ct 3 'P 0 - 

(48) 

#1 , 

^3 ( r ) “ P0 a 3 Tt ~ V 'Pi ’ 

(49) 

^2 , 

?3 " Po a 3 Tt ~ V 'Pi ’ 

(50) 

As before, and \p 2 are chosen so that v 3 and £ 3 vanish and £ 3 is a constant, c 3 . 

expanded in a fashion similar to that used for Equation 40, 

Now 4>i is 

(r)P n m (0) exp im(<fr-qjt) . 

(51) 


nm 


Ifm = othe terms have no time dependence, Equation 49 shows they are harmonic, and H ncl ( r ) is 
proportional to r" (n+1) . Km/ 0,the solution proceeds as for region II. By defining 


,(P) 


P H nm (P) . 
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and a complex magnetic Reynolds number 2a 3c for region III, 


- 2im a. 


+ i a„ 


it is found that 


Z n'm " Z n m [ n ( n + 1 )'°" 2+a 3c ] “ ( 52 ) 

which has the same form as Equation 43. This leads to 

'Pl = j^ P n (0) eXp im (^ - Wt } [ A nml h nnl <P> + A ™2 h nm2 O)] • (53) 

where 

h nmi(-°) = P" 1 £ n expa 3c p , 

h nm2 (P) = ( -1 )" P~ l K eX P ~ a 3c P ■ 

There is a choice of sign for the real part of a 3c . By choosing a 3r as positive, (p) contains 
exp a 3r p which causes it to diverge. Therefore, all A nml vanish and Equation 53 becomes 





^nm2 


(P) P n m ( 0 ) exp im 


(54) 


In a similar manner 


^2 



h nm2 


(P)P n m ( 0 ) exp im (0 “ ojt ) 


(55) 


Contributions by a fixed outside source may be included in Equation 55. Usually Equations 11-13 
then no longer hold and the calculation is somewhat different. 


11 



CONTINUITY CONDITIONS 


On the boundaries the components of B, the tangential component of E, and the normal com- 
ponent of curl B are continuous. The operator A 2 is defined by (Reference 10) 


A 2 


id d l a 2 

inrs de Sln0 de + S in*edtf 


(56) 


so that 


V 2 


1 a 2 r 
r ar 2 



(57) 


and 


[a 2 + n(n + 1)] P n m (^) exp im<^ - 0 . 


(58) 


From Equations 18 and 20 


rB r - - A 2 p 2 , (59) 

and the fact that all boundaries are spherical implies the continuity of all expansion terms of 0 2 , 
as expressed in Equations 34, 47, and 55. In a similar way 

r (curl B) r = - A 2 \p x 


implies the continuity of the terms of \p x . On a spherical boundary curl n p v is also continuous. 

By Equation 18 so is curl curl x\p 2 and it may be shown from this that d\p 2 Jdr is also continuous 
across the boundaries. Finally, by using the above results with Equations 21 and 23 P 0 is also 
found to be continuous. Thus four conditions are to be met by each of the two boundaries and by 
eight sets of undetermined coefficients. Because the source of the field is poloidal, the continuity 
of the components of p 2 and dip 2 /dr makes it possible to evaluate \p 2 independently. Next \p t is 
found, by using Equations 30, 37, and 48 to express the continuity of \p 0 . Finally P Q is obtained by 
using the same three equations. 

As cr 1 and < j 2 increase without limit, the solutions tend to approach those obtained before by 
assuming E = - (v x B) . In general, the solutions diverge in region III unless a 3 is kept finite. 

This may be interpreted as meaning that, in a system having stationary and rotating ideal con- 
ductors in contact, infinite currents will be excited by unipolar induction. If the hydromagnetic 
equation were used this divergence would not occur, because infinite forces would act on the flow 
and distort its pattern. 
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In the example that follows this divergence is avoided in a different way, namely by the 
assumption that in region III conductivity is arbitrary (finite or infinite) along lines of force, and 
zero across them. 

A WORKED EXAMPLE 

As a simple illustration of the preceding discussion, the case will now be considered for a 
field source of dipole of moment M, located at the origin and aligned with the axis of rotation. This 
source produces an axisymmetric field and therefore the variables <f> and t as well as the expan- 
sion index m are absent. It is found convenient to use dimensionless units and to define 


r 


p 

1 o 
ii 


^0 

^01 

u R 0 ’ 

2a x 

= (i 0 a lU R 0 


The scalars defining the field then obey the following equations. In region I 


V 2 ^ = 0 , 

(31) 

V 2 ^ = 0 , 

(32) 

and from Equation 30 


CO 2 1 $ 

'/'oi = - ^ P sin 0-^0- - ^ + c n • 

(60) 

In region II 


<3 

to 

1 

*1!? 

•§1* 

■*> 

u 

o 

(61) 

__ 0 2a d 

V ^2 " — P) ~ 0 - 

(62) 

*" V'oi = 'PlP - 'E'ipty'l p) + C 21 • 

(63) 

And in region HI, if an isotropic conductivity <r 3 is assumed, 


V 2 -/-! = 0 , 

(64) 
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(65) 




V 2 i /; 2 = 0 , 


^01 



( 66 ) 


As mentioned at the end of the previous section, in order to prevent divergence as cr 3 -co and also 
to take into account, at least partially, the anisotropy of conductivity, it will be assumed that in 
region III the conductivity is o- 3 for flow along lines of force and zero for flow across them. Then, 
in region III 


B X curl B - 
B * curl B = 

These equations are nonlinear; however, it will be shown that a solution having the required form 
near the boundary is obtained by taking 


0 , 

M 0 ^3 (E • B) . 


(67) 

( 68 ) 


from which 


curl B - 0 , 


(69) 


E • B = 0 . 


(70) 


By using Equations 18 and 20 and applying the treatment of Equation 24, Equation 69 gives 


v 2 ^ 2 = 0 , 
*, = 0 . 


The unknown scalars are expanded: 

For region I, as in Equations 33 and 34, 


(71) 

(72) 


'Pi = 2Z a "^ P " (0) ’ 

'Pi = P" P n + PP ~ 2 COS 9 ■ 

For region II 

'Pl = + 3 " 2 g n2^)] P n( g ) > 


(73) 


(74) 


(75) 
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HI 


*2 = 22 [ b -> gni (p) + b " 2 ^ (p) ] p " (e) 


And for region III, with p 1 = R^Rq , 


= 2 ^ b .(t-) p .< 9 >- 


(76) 


(77) 


To begin with, 0 2 is derived. For n / 1 the continuity equations of </> 2 and dp 2 /dr yield a set of four 
equations with unknowns b n , b nl , b n2 and B n . Because of the form of the source term in Equa- 
tion 74 these equations are homogeneous and all coefficients vanish. For n = 1, by Equation 44, 


g„ (P) - ( a P~ 1 “ P~ 2 ) exp 2ap , 

g]2 (P) = a P _1 + P -2 ■ 

Denoting differentiation with respect to p by primes gives 

b i + P ~ b n 6n ( 1) + b 12 g u (l) , 
b , - P = b u gjj ( 1) + b 12 g 12 ( 1 ) , 






(78) 

(79) 


(80) 


B 1 - bllBl^PO +b 12gl2(Pl) ' 

-2Bi = b u Pl g,', (p x ) + b 12 Pj g 12 (Pi) • 

J 

The complete solution is rather lengthy. If a » l, which represents the "high conductivity" case, 
and p l » 1, in region I 


+2 


J3 cos 6 


(4 p + p" 2 ) • 


(81) 


in region II 


and in region HI 


(if ) cos e [l - (2apj) -1 exp - 2a (pj - p )] 


(&Pl , 

*a * • cos 9 ■ 


(82) 


( 83 ) 
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The poloidal field in region HI is thus that of an axial dipole, whereas in region n 


) cos 6 p~ 2 , 


(m\ 

3 = R < 

r \ K ° / 

= ('2'kf/Oj) sin 0 exp - 2a ( Pl - p) , 


which may be compared with Equations 6-8. In region I 

B r = cos0(2^ 3 + l) , 

B e = (^) sind {p ~ 3 ~ X ) ■ 

The lines of force of this solution are given in Figure 2, for the limit of high conductivity. 



Figure 2— Lines of force for the poloidal component, in the limit of very high conductivity. 


Because of Equation 72, 4 \ is expressed by only three sets of undetermined coefficients and 
the fact that \p 0 is continuous when r = R x does not have to be invoked. By Equation 30, in region I 


( b i p2+ -p)- 


e ~ 


Id. . 


+ c, 


( 86 ) 


Because 


sin 2 9 3 T( p 0 - p 2) - (87) 

only for n = 0,2 do the coefficients of satisfy inhomogeneous equations and therefore differ 
from zero. The monopole terms do not contribute to the magnetic field, so that only the quadru- 
pole part has to be evaluated. In region II 


from which 


where 


g 2 ! (p) - (a 2 p" 1 - 3ap~ 2 + 3p 3 ) exp 2a p 

g 22 (P) = a. 2 p~ l + 3ap“ 2 + 3p~ 3 , 


a 2 3 21 ®21 * 3 22 ^22^) » 

® ” 3 21 g 21 (Pi) + 3 22 g 22 (^l) 
2 a x 7 - 3 a 2 = a 21 + a 22 k 2 , 


2 a 


( 88 ) 

(89) 


(90) 


A_l = - 3 a j ^1 - exp 2 a , 

A- 2 - 2 aj ^a 2 + 3a + "2 + 3a -1 ^ . 

This again tends to produce cumbersome expressions. When a » 1 and p x » 1 , in region II 

xp x % yp~ x [l - exp 2a (p - j P 2 (0) + f(r) , ( 91 ) 

and in region I 


i// 1 % 7/0 P 2 (< 9 ) + const 


(92) 
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In region n, by Equation 19, 


% 3ryp 1 £l - exp 2a (p - sin 0 cos 6 , (93) 

which should be compared with Equation 7. 

Finally, \p 0 is derived. Equation 70 implies that in region HI the lines of force of B lie on 
equipotentials of \p 0 . B in that region was found to be a dipole field and its lines of force are 


sin 2 8 


cons t . 


Thus in region III 


<A 0 ( r > e ) - <A 0 


sin 2 6 


sin 2 


(94) 


By Equations 44 and 45, or by direct calculation of a 0 -independent solution of Equation 61, the 
monopole term of ^ is, in region II, 


\jj x (mono) - a 01 yO 1 exp 2a p + a Q2 p 


which with Equation 63 gives 


<A 0 1 ( mono ) a 02 + ^21 

By assuming as before that a» 1 and p x » 1, Equations 63 and 91 give, in region II, 

' A 01 “ (^) + ^*22 * 

Combining this with Equations 87 and 94 shows that 

C 22 “ y • 


Thus in region II 


’Aoi 


- ~2 y sin 2 8 , 


and in region HI 


’Aoi 


~ ~2 y sin 2 6 



(95) 


(96) 
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It should be noted that in region II fi 01 does not depend on r. The equipotentials thus tend to be 
(with large a and p x ) cones of constant 9. 

The potential in region I is determined in the same manner. If a i » i > the term containing 
</>! in Equation 86 may be neglected and 

'/'Oi = (ir)(§ P 2 +/ 5 ' 1 ) sin2 ^ • (97) 

The approximation involved here is the same as for Equation 2, so that Equation 70 holds in all 
regions and the magnetic lines of force lie on equipotential surfaces. Thus Figure 2 may also be 
viewed as a cross section of the equipotentials of the electric field. 

It is interesting to note that, in the limiting case of high conductivity, the electric potential 
in region II, especially near the equatorial plane, is different from that at infinity. In particular, 
the equatorial plane is an equipotential in which, by Equation 97, 


0o 


2 wR o P 


(98) 


The quantity /3/r 0 , which is of the order of the field at r = R 0 , will be taken as 1 gauss; Equa- 
tion 98 then gives a potential of about 2 x 10 8 volts. This result may have some connection with 
the modulation of cosmic radiation by the solar activity cycle which resembles that produced by 
an electric field (Reference 14); however, the value of deduced here is too small by a factor 
5-10, and it should be borne in mind that when the solar dipole reverses its direction, as observed 
in 1958 (Reference 15), <// 0 reverses its sign. 

To the preceding example may be added a homogeneous "interstellar magnetic field" B 0 
which, in order to preserve symmetry, will be assumed to be parallel to the rotation axis (for an 
arbitrarily directed B 0 the calculation is more involved). Such a field can be represented by a 
poloidal potential 


«A 2 - ~2 B 0 r cos ^ 


(99) 


so that Equation 77 is replaced by 


02 



P n <*) 


2 B o R o ^ cos e 


( 100 ) 


The inclusion of B 0 causes a term RjBq/ 2 to be added on the left of the last two of Equa- 
tions 80. When these equations are solved, it turns out that b 12 (the important coefficient in 
region II) is modified by a factor of 


1 


+ 



exp 2a ( 1 - /Oj) , 
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which generally will differ very slightly from unity for large a and p x . Thus </> 2 , and consequently 
\p l and ^ 0 , are only negligibly affected in region I and in most of region II. Note that it is quite 
possible for the outlying interplanetary field to be much weaker than the surrounding interstellar 
one. The solar wind then scoops out a cavity in the interstellar field, as first suggested by 
Davis (Reference 16). 
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